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I. INTRODUCTION 

The Berezinskii-Kosterlitz-Thouless (BKT) transition is well-knovifn as an infinite-order phase transition Q. The 
correlation length ^ has an essential singularity at the critical coupling parameter gc 

^^eMMg-gm. (1) 

with the critical exponent cr = 1/2 or 1. In c = 1 conformal field theory(CFT) [||, there are infinitely many models 
where infinite-order phase transitions can occur. Any of them shows the same universality as the BKT transition. 

One observes a different from 1/2 or 1 in some c > 1 CFTs |^-^. Recently, a model of a quantum spin chain, whose 
long-distance behavior is described by level 1 SU{N) WZW model, was studied by Itoi and Kato They pointed 
out that an infinite-order phase transition with a critical exponent a = occurs by an SU{N) symmetry-breaking 
marginal operator. In the iV = 3 case, this corresponds to the gapless Haldane gap phase transition in a spin 1 
isotropic antiferromagnet in one dimension. In a problem of dislocation-mediated melting, some curious numbers 
were observed by Young, Nelson and Halperin Q. They obtained ct = 1/2 for a model on a square lattice, a — 2/5 
for a simplified model on a triangular lattice and a non-algebraic number a — 0.36963... for a generalized model. In 
Ref. 1^ , Bulgadaev studied topological phase transitions in c > 1 CFT with non-abelian symmetry, where non-abelian 
vortices play an important role. They belong to special classes of infinite-order phase transitions and several series 
of a dependent on the symmetry of the system, were found. Though there have been several studies for those some 
different types of models with infinite-order phase transitions including BKT type, the universality classification by 
this critical exponent still remains a challenging problem. 

In this paper, we study the universal nature of the critical exponent a in infinite-order phase transitions. We show 
that the critical exponent a is determined from the operator product coefficients of the marginal operators which 
cause the infinite-order phase transition. It is shown that a marginally irrelevant operator can also affect the value of 
the critical exponent a. 

The critical exponent of the correlation length is extracted from a long-distance asymptotic form of running coupling 
constants, whose leading term is determined by the motion of the coupling constants near a fixed point. In an ordinary 
finite-order phase transition, we linearize the renormalization-group equation (RGE) around the fixed point and can 
derive the exponent exactly. Namely, we can show that the inverse of the exponent is equal to the maximal eigenvalue 
of the scaling matrix defined by the derivative of the beta function at the fixed point. One does not have to solve the 
differential equation exactly in order to obtain the exact critical exponents in this case. In the infinite-order phase 
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transition however, the scaling matrix vanishes since the phase transition is driven only by marginal operators. So 
far, one has had to solve the differential equations explicitly to obtain the critical exponent, although RGEs with 
multiple variables are generally non-integrable due to their non linearity except for some fortunate cases like the BKT 
transition. This difficulty is one of the reasons why the universality classification of infinite-order phase transitions 
by the critical exponent a in eq.(^ has never been successfully done. 

In order to resolve this problem, we apply another renormalization group (RG) method developed in Refs. j^J^ to 
studying the long-distance asymptotic behavior of the solution of the original non-integrable RGE. This RG method is 
starting to be recognized as a general tool for asymptotic analysis. Chen, Goldenfeld and Oono Q introduced the idea 
of RG to singular perturbation theory and gave a unified treatment. According to Bricmont, Kupiainen and Lin 
the RG transformation for a partial differential equation is defined as a semi-group transformation on a space of initial 
data, which is generated by a scaling transformation combined with time evolution. Koike, Hara and Adachi used 
this general method practically in the study of the critical phenomenon in the Einstein equation of the gravitational 
collapse with formation of black holes [|| . Tasaki gave a pedagogical example for the RG transformation, where the 
equations of motion in Newtonian gravity were analyzed ||^ . 

In SecH we reuse the RG transformation of Ref . jgj , which enables us to calculate the critical exponent a in cq. (|l|) 
without solving the nonlinear differential equation explicitly. The new RGE regards the straight flow line solving 
the original RG E a s a fixed point, where the derivative of the beta function in the new RGE has non-zero value in 



general. In Sec, III, we show that the inverse of the maximal eigenvalue of the scaling matrix derived from the new 



RGE gives the critical exponent a. In Sec. IV, we also study asymptotic behavior of the running coupling constants 
in a massless phase and extend the well-known formula for a logarithmic finite-size correction to the case of multiple 
running coupling constants. In Sec.^ we exhibit several nontrivial examples motivated by antiferromagnetic quantum 



spin chains. Finally, we give a summary and discussions in Sec. VI. 



II. RGE FOR RGE 
A. Formalism 

Let us begin with the RGE for a given set of n marginal operators 

where g = (51, • • • ,gn) is a set of coupling parameters and t = log/ with / being a length scale parameter. Since the 
operators are all marginal, the right-hand side is expanded as 

H(g) = ^C^^g,:5,+0(g3), (3) 

where C^^ is proportional to the operator product coefficients of the operators. First we neglect the higher-order 
terms 0{g^\ and later we discuss irrelevance of those neglected terms. 

In general, we find several critical surfaces where the RG flow is absorbed into the origin. A phase transition occurs 
if the initial coupling constants cross one of the critical surfaces. These critical surfaces divide the coupling parameter 
space into several regions which are phases. In the next section, we consider one massive phase surrounded by a set 
of critical surfaces, where there are several marginally relevant coupling parameters. In this region, we have a finite 
correlation length, which becomes larger as the coupling parameter approaches the critical surface. 

We are going to study the long-distance asymptotic behavior of solutions for the RGE (||) which is non-integrable 
in general. To this end, let us introduce the RG method explained in the introduction. We define a renormalization 
group transformation on n — 1 dimensional sphere which forms a set of initial values. We denote the solution g of 
eq.(|^) with the initial condition ap — (floij ■ ■ ■ : O'Qn) as 

9(i,ao), (4) 

namely, g(0, Oq) = Qq. The function g{e^t, ap) is a solution of the RGE as well, because of its scale invariance. 
Let S be the n — \ dimensional sphere whose center is at the origin with the radius |ao| = oq. We define a new 
renormalization-group transformation IZr : S ^ S as follows: 

TZrao = p7g(s(T), ao) = a(T). (5) 
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Note that TZt has a semi-group property: 

7^rl+r2 = T^T2 ° "^Ti- (6) 

The meaning of TZt is the fohowing: first, choose t. Then move ao along the solution g(t^ao) during the time s(r). 
Here s(r) is determined by the condition g{s{T),ao)e'^ € S. See Fig.^ 

g{s{T),ao) = a{T) 




s{T),ao) 



FIG. 1. Illustration for TZt and the beta function defined in eq.(^). For simplicity, we take n — 2. The dashed line 
represents the tangent space at a(r) £ S. 

Next let us derive the beta function for TZt- Noting that V{g) is quadratic, we have 

da T^r, , \\ds 

— = a + e V{g{s,ao)) — 

CLT ClT 

ds 

= a + e-^Via)^. (7) 
or 



The length-preserving condition 



a.-=0 (8) 



leads to the following differential equation for s(t): 



dr a ■ V{a) 

with the initial condition s(0) = 0. Inserting cq.(||) into eq.(0), we obtain the beta function for TZt- 



(9) 



Note that f3 can be written as 



^da^^ a.a-V{a)-V.{a)al^ 
dr a - V{a) 



a ■ V(a) 



where P is the n x n matrix that projects V(a(r)) onto the tangent space at a(r) e S: 

P,^ia)^6.,~^. (12) 

For later use, we derive a polar coordinates representation of the new RGE. Employing polar coordinates, a e is 
expressed as 

CTL—l n— 1 TL—1 \ 

ao Y[ sin^Q, ao cos6'i JJ^ sin^Q, ao cos02 sin6'Q,, • • • , oq cos0„_i J . (13) 
a=l a=2 a=3 / 
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Since {da/d9a}i<a<n-i are orthogonal to each other, we can make the basis {ea}a orthonormal on the tangent space 
at a € S* by an appropriate rescaUng 



Ba = fa{a) /a(a) = 



da 



(14) 



Then, 



Pa= 13 ■ ea ^ ^ ■ Ba = faia)^^, (15) 

ar dT 



which leads to the RGE in polar coordinates 



= (/„(a))-^^„(a). (16) 



Returning to the coordinate- free representation eg. (pi]), let us find a fixed point of the new RGE (p^). The nature 
of the new RGE near the fixed point determines the universal behavior of the infinite-order phase transition, as in 
the ordinary RGE for a finite-order one. Near a fixed point, a(T) moves slower as its trajectory tends to a critical 
surface. This implies that the time a(T) spent in a neighborhood of the fixed point is a singular function of the initial 
condition ag. This singularity can occur only at fixed points of the new RGE, which allows us to analyze its singular 
behavior by a linearization near the fixed points. 

From eq.(|ll|), one finds that a is definitely a fixed point if P{a)V{a) = and a ■ V(a) ^ 0. In this case, since 
V{ka) is parallel to a for all real numbers fc, a is on a straight flow line of the original RGE (||). Straight flow lines are 
put into two classes. If an arbitrary point a on a straight flow line satisfies a ■ V(a) < 0, it is said to be an incoming 
straight How line because a is carried toward the origin in time evolution. On the other hand, if a ■ V{a) > for all 
a on a straight flow line, it is called an outgoing straight How hne. If a fixed point a of eg. (p^) is on an incoming 
straight flow line, —a is a fixed point on an outgoing straight flow line. 

What happens if P(a)V(a) = a ■ V{a) = ? In this case, V{a) itself vanishes. It means that a is a fixed point 
of the original RGE (||). Moreover, since V is homogeneous, ka is also a fixed point for all fc S R. Namely, the 
original RGE (||) has a fixed line in this case. If the original RGE (^ has this fixed line, a point on the fixed line has 
a non- vanishing scaling matrix even though the coupling constants arc all marginal at the trivial fixed point g — 0. 
Therefore, we can directly analyze the original RGE near a point on the fixed line and can show that the phase 
transition generally becomes of finite order in this case. 

Here, we give a couple of remarks on the global nature of the new RG transformation TZ^ defined by (^. First, 
there could be a turning point g where V{g) ■ g — with V{g) ^ 0. Let 

log^... (17, 

Although egs.(||) and ( po|) cannot be defined at r = f, it is obvious in a geometric sense that a(f) and s(f) are 
well-defined. For example, in Fig.^ ^■(t^) = and sir) are determined as the definite time g(t^ ag) spent during the 
journey from to g. 




FIG. 2. An example in the case of a flow having a turning point. Here we take n = 3. The gray lines represent solutions of 
the original RGE (^), while black ones on S for the new RGE (0). Here a* (a'*) on the incoming(outgoing) straight flow line 
is a fixed point of the new RGE. 
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Second, if g{t, Oq) has turning points, a(r) and s{t) become multi-valued with respect to r. For example, in Fig||, 
g(t,ao) has a turning point at g. Suppose that = Ig^^^l = b and choose t = log a/6 = tq. Then Ttro'^o has 

two images, e^^gr^^^ and e^°g^'^\ In this case we distinguish the images as a(^)(ro) and a^'^\To). Similarly, s{tq) also 
has the same multiplicity, which is distinguished in a similar way. In Fig.^, the image of TZt starting at ao reaches 
a branch point e^g. We denote the solution from ag to e^g by a'^^^(T). The remaining part is called a(2)(T). Both 
a^^-'(T) and a'^-'(T) are absorbed into the branch point ge^ , which indicates the fact that g gives a minimum distance 
from the origin. If a turning point corresponds to a maximum distance, the two solutions will escape from the branch 
point. 



B. Example — the 2D classical XY model 



Here we exhibit the new RGE for the 2D classical XY model as an illustrative example. The original RGE for the 
XY model is given by [|ol 



dgi 
dt 

dg2 
dt 



V{g) 



~9l 
'9i92 



(18) 



with g2 > 0. Let us first look at the phase structure from eg. (|l8|) . See Fig.||(a). 

(a) 



92 
A 



92 

i 



(b) 




<^9i 

FIG. 3. (a) Flow generated by the original RGE of the XY model, (b) Flow of the new RGE ^ with eq.(|l|). Fixed 
points are denoted by the black circles. The new RGE is not defined at the white circles which correspond to a fixed line or 
turning points in Fig.P(a). 

The RGE (|l^ has two straight flow lines gi±g2 = and one fixed line 172 = 0. It is well known that each point on 
the fixed line corresponds to the theory of a 2D massless free boson that is parametrized by a compactification radius 
of the boson field [g). The shaded region in FigJ^(a) is a massless phase since flow in this region is flnally absorbed 
into a point on the 52 = fixed line. The incoming straight flow line, gi — 52 = with (72 > forms the phase 
boundary. As an initial coupling approaches the phase boundary from the massive phase, the correlation length ^ 
becomes divergent.^ 

Now we turn to the new RGE for the XY model, which is given by eqs.(p^ and ( [l^ ) with the condition al + a^ = a§ 
{0-2 > 0). It is explicitly represented as 



dai 



da2 



4) 

2k ("2 - «?) 



in cartesian coordinates. Alternatively, using polar coordinates a 
becomes 



(19) 



(aosin^, oocos^) (-7r/2 < 6* < 7r/2), the RGE 



^ The correlation length also diverges when the initial coupling constants tend to a fixed point on 172 = with gi < 0. However, 
the scaling matrix at this point does not vanish and the ordinary finite-order phase transition takes place. Since our interest 
is focused on an infinite-order phase transition, we do not consider that case here. 
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— = -cot 261 (20) 



owing to eq.(|l6|). 

Next we find fixed points of the new RGE. Solving 



P{a)V{a) = 




we have a — (ioo/v^, ao/\/2), (±ao, 0). Evaluating V{a) ■ a at those points, it turns out that (oq/ V2, ao/ V2) is 
on an incoming straight flow line while (ao/\/2, — ao/V2) on an outgoing straight flow line. The remaining points 
(±ao,0) are found to be on a fixed line. Note that V{a) ^ and V{a) ■ a = when a = (0, ao). This means that 
(0, oo) corresponds to turning points on trajectories generated by the original RGE. The flow of the new RGE cannot 
be defined at the points (0, oq) and (±ao, 0). Since the example is a simple two-parameter system, we can understand 
qualitative aspects of the global flow in the new RGE. See Fig.|(b). 

In the last part of the next section, we will continue the analysis of this model and derive a from the beta function 
in cq.(|l^). Before performing that, we need to know a representation of the correlation length in terms of the new 
RGE. 



III. CRITICAL EXPONENT OF THE CORRELATION LENGTH IN A MASSIVE PHASE 



In this section, we explain how to evaluate the critical exponent a of the correlation length in a massive phase from 
the beta function (10). 

We first define a correlation length ^(oo) by the following formula: 

|g(logf(ao),ao)| = 1. (22) 

Namely, log^(ao) is the time g(t, Oq) spent in the perturbative region. We note that ■^(ao) defined above changes as 

e*^(g(t,ao))=e(ao) (23) 

under the original RG transformation, which should be satisfied by an intrinsic length scale of the model. The 
differential form of this equation is obtained by eq.(^ 

which is well-known as an equation for an invariant length scale £,{g). Further, eq.(]2^) is a natural generalization of 
the correlation length used in the 2D classical XY model. 

We consider the case where the running coupling constants are in the perturbative region \g(t,ao)\ < 1. In this 
section, we study in particular the long-time asymptotics of a flow that once approaches the origin and then leaves for 
the non-perturbative region \g(t, ao)| > 1, as the flow in the region g2 > \gi\ in the XY model. Generally, a quadratic 
differential equation such as eq.(H) admits a flow qualitatively different from that investigated here. In section VI, we 
discuss such an exceptional case. 

Next, let us represent £,{ao) by the solution of the new RGE(|l^). Dcflne tr by e-''^\aa\ ^ 1. From the definition 
of s(t), we obtain 

loge(ao) = s(TR) = dT — . (24) 

Using the differential equation (^) on the right-hand side, we obtain the integral representation for ^(oto). Since a 
flow treated in this section has a turning point as shown in Fig.|^, the correlation length is represented via 

'^^^^-^^ = -1 ^^a(i)(r) .V(a(^)(r)) " /, ^ '^a(^)(.) .V(a(^)(r)) ' ^''^ 
Employing the integral representation, we argue that the leading term of ^ is given by 

loge(ao)-e^ (26) 
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if f in eq.(^^ is sufficiently large. Even though the integral near the turning point seems to diverge, it is only 
apparent as discussed in the previous section. The first term in the right-hand side of eq. diverges due to in 
the integrand when f goes to infinity. The second term contributes to the correlation length with the same order as 
the first term. Hence, we can evaluate C(ao) by eq.(^6|), which translates singular behavior of ^(ao) into that of f. 

Next, we evaluate the divergent f by using the polar-coordinate expression eq.(]T^) of the new RGE. ft is obvious 
that f diverges if the initial coupling constant ag is on an incoming straight flow line. It implies that f grows when 
a(r) passes near a fixed point a* on the incoming straight flow line. 

Suppose that a(r) goes through a neighborhood U oi a* , as shown in Fig.|[ 




critical surface 



FIG. 4. Flow in the new RGE ( |10[ ) near the fixed point a* on an incoming straight flow line. Dashed line represents 
one-parameter family of the initial value ao{T). 

The scaling matrix of the beta function ( [l^ ) does not vanish at the fixed point in general and then the beta function 
can be linearized in U . That is, 



7=1 
n-1 

7=1 



(28) 



where S9^ = 9^ — 9* with {9*}^ representing the fixed point a* . The scaling matrix Aay{a*) describes the large f 
behavior because a(T) spends a long time in U. If the scaling matrix is diagonalized with eigenvalues ba by a new 
coordinate {9'a}a, the new RGE becomes 



dT 



t - b^d9' 



(29) 



The solution is 



t{r) = 59'Me'"'^^-^'^. 



(30) 



^ In fact, there exists a component of V which does not vanish at the turning point, say, Vi. We parametrize the flow of the 
new RGE by ai instead of r near the turning point and change the integration variable in eq. (|2^). The measure changes as 

dr ^ doi 

a{r)-V(a{T)) aia ■ V (a) - a^Viia) ^ ' 

from eq. (^o|). The denominator of the right-hand side does not vanish near the turning point. Therefore the denominators 
in the integrand in eq.(|2^) does not contribute to the divergence of (,{ao) and can be replaced by a certain constant when we 
evaluate the leading divergence of the integration in eq. (pSj) . 
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We take U as the (n — l)-dimensional cubic box, whose side is 2e and center is a* . If the scaling matrix has the 
unique relevant mode O'^, a{T) spends time 



7^ log 

Ol 



(31) 



<50i(^o) 

in U . Here we have supposed that a{T) reaches U aX t — tq. Now we vary the initial value ao by one parameter T 
and assume that ao{T) intersects a critical surface transversally at T = Tc- See Fig.^. 

As the initial value ao{T) tends to the critical surface, |(56''j^(to)| gets small. It implies that 66[{to) is expanded as 

5e[ {to) = const. (T - T,} + O ((T - Tj') . (32) 

Since a(T) spends a finite time outside of U, the divergent behavior of ^(ao) is determined by eqs. ( ^l| ) and (|3^). 
Thus we get 

logaao)^e' ^IT-T^^'^ (33) 

which means 

„1. ,34) 

This quantity docs not depend strongly on the choice of ao(r), so that cr is a universal quantity in this sense. 
It is quite useful to find a relationship between the scaling matrix A{a*) in eq.(|2^) and the n x n matrix 

S,,(a*)^g(a*) (35) 

for practically computing the eigenvalues {&q}q. In appendix we will show that 

A{B) = A(yl) U {0}, (36) 

where A(M) is the set of eigenvalues of a matrix M. It should be noted that B{a*) = B{—a*) since /3(a) is an odd 
function. This means that the scaling matrix A{a*) has the same eigenvalues as A{-~a*). 

Now we deal with the 2D classical XY model again and show how to derive a by our method. The original RGE is 
given by eq. (|l^) . We saw in the previous section that a* = (oq / aq / V2) is a fixed point on the incoming straight 
flow line. The matrix B defined by eq.(|35|) is 

i.e., the scaling matrix A{a*) has the eigenvalue 2, according to eq. (^6|) . Alternatively, A{a*) is directly computed in 
this case. From eq.(pO|) we get 



Aia*) = 



cot 261 = 2, (38) 

ei=7r/4 



as expected. 

Following our result cq.(^4|), we get 

a = 1/2, (39) 

which is well-known as the BKT universality. Originally it is obtained integrating the nonlinear RGE ( p^ ) explicitly 
1^^ . In contrast, according to our approach, we can reach the same result in an algebraic way. As we will see in 
section^ it can provide a even in the case when an original RGE(^ is not integrable. 

In the above example, the fixed point a* has a unique relevant mode, so that we can apply the result (|3^). If there 
are multiple relevant modes in the scaling matrix, we can observe other relevant exponents l/t'2, 1/^3, • • ■ > in an 
appropriate fine tuning of the initial parameters. 

Finally, we discuss the irrelevance of the higher-order terms in the original RGE (||). Here, we assume that we 
acquire no extra fixed points by taking into account higher-order terms. If we have higher-order terms, the RG 
transformed coupling with r obeys a different equation because of their scale breaking nature. The scaled coupling 
g'it) = e''g{e'^t,ao) obeys 

^ = W) + 0(e-V'). (40) 

Note that the higher-order term becomes smaller and the RGE takes the scale invariant form asymptotically. Therefore 
higher-order terms are irrelevant to determine the critical exponent. 
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IV. LOGARITHMIC DEPENDENCE OF MULTIPLE MARGINALLY IRRELEVANT COUPLING 

CONSTANTS 



So far, we have studied solutions of RGE (g) in a massive phase. Our method is also applicable to studying 
asymptotic behavior of a solution in a massless phase. In this section, we study the logarithmic dependence of the 
multiple running coupling constants in a massless phase. 

It is important to clarify finite-size corrections in a system with marginally irrelevant operators. For example, a 
numerical simulation in a spin system can calculate energy levels only for small degrees of freedom. A theoretical 
formula for the finite-size correction is useful to extrapolate numerical data to those in the infinite system. If the 
system can be described in a critical theory with marginally irrelevant perturbations, physical quantities acquire 
logarithmic finite-size corrections. Here we are interested in a system with a finite volume described by a theory 
with marginally irrelevant coupling constants g, where D is the space dimension. Consider the situation where we have 
a critical theory with g = Q which describes the system with an infinite volume. In the system with a finite volume 
, we can calculate physical quantities with a finite size correction in the critical theory with a small perturbation 
of the coupling g obeying RGE (|^). If we have an initial coupling ao at a lattice spacing 1, the running coupling at 
the scale L becomes g(logL, Og), where g(oo, Oq) = 0. In the case of a single marginally irrelevant coupling constant 
5, the original beta function is given by 

V{9) = Cg^ + 0{9^), 

where C is a universal constant in the sense that it is independent of an initial value. The running coupling constant 
with an initial condition ap has the following solution 

rg{logL,ao) 1 ]_ 

~Jao nff) ^ Cg{logL,ao) ~ C^o' ^ ^' 

In this solution, we have a well-known universal expression 

.(logL, ao)^^^ + o( -i-, 'f^) . (42) 
Clogi V(logi)^ (logi)V 

The leading term is independent of the initial coupling ao, and therefore this formula is useful to fit numerical or 
experimental data of the system with a finite size. For example, in one dimensional quantum spin systems with 
marginally irrelevant perturbations, logarithmic finite-size corrections to the ground state energy 

/\Eo = - — \c + —7 + O 



QL\' (logL)3 V (logi)'* ' (logL)4 



are calculated from this formula (^) ]11| , |12| , where c is the central charge and A is determined from C. Since c and 
A are universal constants, we can compare c and A to numerical (experimental) data and obtain a clue whether a 
field theory that derives eq.(^) is truly effective or not. Therefore it is important to derive a formula corresponding 
to eq.(^2|) where there are multiple marginally irrelevant couplings. In this section, we shall show that this universal 
nature holds in this C&rSC, clS well. 

As we mentioned above, we examine the case where the all coupling constants are marginally irrelevant, so that 
fiow of eq.(^ is absorbed into the origin. In this case, there are no turning points on the flow, which implies that 
the transformation TZtCio defined in eq.(^ is single- valued with respect to r. Therefore we can write down a formula 
similar to eq.(|25|) as 

logi = - rdr' (43) 
Jo air') ■ V{a{T')) 

The running coupling constant g{log L, ao) is obtained by 

Sf(logL, ao) e-^a(r). (44) 

In order to derive the logarithmic dependence of ^(logL, ao), we first solve eq.(^3|) for t when L is sufficiently large. 



Then we apply the result to eq.(44). 

As we have seen in the previous section, when we take L sufficiently large, the contribution from a neighborhood U 
of a fixed point a* on an incoming straight flow line dominates in the integration of eq. mq) , which can be evaluated 



from the hnearized new RGE in U. Suppose that a(T) enters J7 at t = tq. Eq.(M3) becomes 
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ale7' 



for large L. 
Writing 



logi^-/ dr' \ (45) 



a(T) = a* ^ba{T), (46) 



(5a(r) in the polar-coordinate representation obeys the linearized RGE ( p9|) in U . Its solution has the following 
asymptotic form for large r 

'^"W - E ^^^^W - ^5^'i(ro)e^^'---\ (47) 

Q — 1 ^ 1 

where &i < is the maximal eigenvalue of the scaling matrix A{a*) defined in eq.(^8|). We expand the integrand in 
eq.(E|) as 



,2 „r' „2 ^t' 

and calculate the right hand side of (^). First we compute the leading-order contribution. The leading integration 
is easily performed as follows: 

\ogL^~ Tdr' f^f. . =- ^ g. . (e--e-»)c.- ^ g e^. (49) 
a* ■ V[a*) a* ■ V(a*) a* ■ V{a*) 

Since a ■ V{a), which is a cubic function of {a^}, is negative at a*, we can write 

a* ■ V{a*) = -Cal, (50) 



where C is a positive constant defined by 



C = -e* - Vie*) (51) 



with e* = a* /uq. From eqs. (H), (49) and (p3), we get, in the leading order, 



g(logL,ao) = e-W~^^e*. (52) 

Since e* and C are completely determined by the explicit form of V , the result in the leading order is universal. 

Next, let us go to the next-to-leading term. After evaluating the next-to-leading term in the integral (^5|) with the 
help of eqs.(p7|) and (|4^), we represent e"^ in 1/logL expansion. The calculation is easily performed and finally we 
obtain 

g(logL,ao) = { + ^^i^f^^^ (6i = -1) , (53) 

1 „* I B" 



+ Trae* (^^1 < -1) 



ClogL^ ' (logL) 

where the constants B, B' and B" depend on the initial condition, in contrast to the leading term. The result implies 
that, if — 1 < 6i < 0, we have to take into account the non-universal nature of the subleading correction even though 
0{g^) corrections in the original renormalization-group equations give universal coefficients to this subleading term. 

V. EXAMPLES 

Here, we consider the level 1 SU{N) Wess-Zumino-Witten (WZW) model in two dimensions as a critical theory 
This model has traceless chiral currents J°^(z) and J°'''{z) {a,b — 1, ■ ■ ■ , N) satisfying the following operator product 
expansion (OPE): 
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r\z)r^{o) - 0, (54) 

where ~ stands for equality up to regular terms. Using these currents, we can construct {N'^ — 1)-^ marginal operators 
J°^(z) J'^''(z). In this section, we study models perturbed by some of the marginal operators which are inspired by a 
quantum spin chain • 



A. two-parameter system 

First, we consider a simple two-parameter system which includes the BKT universality as the special case N = 2. 
We define the SU{N) symmetric marginal operator (j)^{z,z) and the symmetry breaking one 0^(2;, z) 

N N N N 

a=l b=l a=l b=l 

which satisfy the following closed OPE 

—2N 

0i(z,z)</)i(O,O)~ ^01(0,0) 
0i(z,z)02(o,O) ^ ^ (0^(0,0) - <^^0,0)) 

ON 

02(z,z)</>2(O,O)~— (/)2(0,0) (56) 
by eq.(|5^). The action integral A of the perturbed theory is 

Awzw + "^01 / — 0*(z,z). (57) 
The OPE formula eq.(|56|) yields the following two-parameter RGE: 

^ = Vi{g)^g,{Ng, + 2g2) 

^ = V2{g) = -52(251 +Ng2). (58) 

In the case of A'' = 2, the RGE reduces to the same form as the RGE of the XY model with an appropriate linear 
transformation, which was extensively studied in sections ^ and III. Here we restrict ourselves to the case of iV > 3. 
The beta function in the new RGE (|l^) for eq.(|8|) is 



/3i(a) = ai 



ai{Nai + 2a2)al 



al{Nai + 2a2) - a^(2ai + Na2) 



a - -a2{2ai + Na2)al . . 

/32(a) ^a2- ^ 2a2) - al{2a, + Na2) ' ^^^^ 

Solving /3(a) = 0, we have the following six fixed points: ±(ao,0), ±(0,ao), and ±(ao/V2, — ao/\/2). Evaluating 
a ■ V{a) at those points, we find that there arc the three fixed points on incoming straight flow lines, (— ao,0) = 
ci, (0,ao) = C2, and (— ao/\/2, ao/V2) = C3. The matrix B{a) in eq.(^5|) at those points becomes 

B(co.(;4,),«,c„.(f S),.,C3,.|±^(;;).(;4), 

which means Ci and C2 are unstable fixed points while C3 is stable for all iV > 3. 
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Divergence of the correlation length is governed by the unstable fixed points and, according to the formula eg. (p4| 

N 



N- 



(61) 



which is identical to that obtained by the explicit solution of the differential equation in Ref. . Note that the result 
eq.(|6^) is also valid in the case of TV = 2 although the scaling matrix cannot be defined at C3: C3 corresponds to a 
fixed line of the original RGE (||) if TV = 2. 

Any theory of c = 1 CFT with marginal perturbations has the critical exponent cr = 1/2 or cr = 1, as is well-known 
| p5[ . This is because the level 1 SU{2) WZW theory is the maximally symmetric theory in c = 1 CFT and because it 
gives the most general theory with marginal perturbation in c = 1 CFT |l^ . The most general theory with marginal 
perturbations describes a quantum XYZ chain with spin 1/2. The infinite-order phase transition occurs at a line of 
the XXX chain. This corresponds to = 2 in our analysis. In the case of c > 1 CFT with marginal perturbation, 
however, we show some new universality classes with non-trivial critical exponents a eq.(|6l|) for N > 2. For example, 
the transition in the case of iV = 3 describes the gapless Haldane gap phase transition with the exponent cr = 3/5 
from the SU{3) symmetric line g2 = in an isotropic spin 1 chain Q. 

The result eq. ( |60|) is also useful when we figure out the qualitative picture of the flow in the original RGE (^8|) . For 
this purpose, we need to know a branch point, which corresponds to a turning point on a flow in eq.(p8|), by solving 
a ■ V{a) — 0. The solution is ±(ao/\/2, ao/V^) for all > 3. A flow in the new RGE changes its direction at these 
points, as we depicted in Fig.||. Combining the result eq.(pO[) and the fact that the scaling matrix at — (i = 1,2,3) 
has the same eigenvalues as at c^, we get the global flow of the new RGE, as in FigJ|(a). It derives qualitative features 
of the RG flow in eq.(58), which are drawn in gray curves in Fig.H(b). We notice that the region gi < 0, (72 > is a 
massless phase, where solutions in the original RGE are all absorbed into the origin along the incoming straight flow 
line gi + 52 = (gi < 0). The incoming straight flow lines passing Ci or C2 form the phase boundary. 





FIG. 5. (a) Flow in the new RGE. The black circles stand for fixed points while the white ones branch points corresponding 
to turning points on the flow in eq.(M). (b) Illustration of the flow in eq.(Ba) derived by the new RGE, which is drawn in gray 



Next, we discuss logarithmic dependence of the running coupling constants in the massless phase. Let us introduce 
new variables {X, Y) = {gi — g2, —gi — (72)- The original RGE has the incoming straight flow line Y ^ {X < 
0) on which C3 is situated. According to our result eq.(|5|), the running coupling constant A'(logi,ao) has the 
leading logarithmic dependence 1/logL, whose coefficient is universal. In contrast, Y{\ogL,ao) has the dependence 
(logL)-i+'' with a non-universal coefficient, where b = {N + 2)/(2 - A^) < -1 for A^ > 3. Hence the l/(logL)2 
contribution that belongs the next-to-leading term in A'(logL,ao) gives sub-leading contribution. This implies that, 
if we can determine g'^ terms in the original RGE (p8|), universal loglogL/(logL)^ dependences can be obtained in this 
example. We remark that the logarithmic dependence of y(logL,ao) is consistent with the result from the explicit 
solution ||]. 



B. three-parameter system 

Here, we consider a non-trivial three-parameter system, an S'?7(2)-invariant marginal perturbation of the level 1 
SU{A) WZW model whose RGE becomes nonintegrable. This model may describe an 5* = 3/2 quantum spin chain 
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around the S'C/(4) symmetric Uimin-Lai-Sutherland model Jl3| , p^ with some SU{2) invariant perturbation. The 
SU{2) transformation is generated by 

Tr J dzJ{z)L + Tr J dzJ{z)L, (62) 

where L = {L^,L^,L^) is the spin matrix in the spin-3/2 representation. Marginal operators invariant under the 
SU{2) transformation is constructed as follows: 



'(z, z) = Tr J2 {-lY''J{z)T,^mJ{z)Tj,~rn, J - 0, 1, 2, 3, 



(63) 



where Tj^m satisfies [L'^,Tj^„i] — j{j + 1)?^,™ and [L^, Tj^m] = niTj^m- Using the tracelessness property of the currents, 
we get 



^0^"(z,z)=O, 

which indicates that there arc three independent marginal operators in (jp , 

cPz 



2 

3=0 



27r 



■(/)^(z,z). 



Employing the OPE ( p4[ ) and the normalization condition Tr^TjmT^ 
product expansions 



j'.m' 



(64) 

• • , 0^ . Here we consider the perturbation 

(65) 

5jj'5mm', we find the following operator 



cjP(z,z)<t>\m^T1^4>\m 



0i(z,z)02(O,O) 

02(z,z)02(O,O) 



/(z, z)0i(O, 0) ~ -:-(-|0O(O, 0) - 0) - ^(/.2(0, 0)) 



1 , 1 
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<^°(z, z)0^(0, 0) ^ —(--00(0, 0) + -4>^{Q, 0)) 
<^i(z, z-)0i(0, 0) ^ ^(-1^0(0, 0) + ^01(0, 0) + ^<^2(0, 0)) 
i,(-^0O(o, 0) - i0i(0, 0) - 0)) 



10' 



r(30O(O,O) + 302(0,0)). 



(66) 
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The RGE derived from the OPE is 

dt 
dgi 
dt 
dg2 
dt 

We find that the new RGE for eq.(|67|) have the following seven fixed points on incoming straight flow lines: 

T 



Voig) = 


-5o' + 


3 .go .gi 






5 


25 


Vi{g) = 


5o5i 




5i 52 


2 


10 


2 


V2{g) - 


3 .go gi 
10 


3.gi2 
10 


50 52 
2 



(67) 



ci = (ao, 0, 0), C2 = I 0, ) , C3 



C5 



,V2' ' V2 
'3ao 5ao 2ao / 2 



3ao ^ Co 



10 



/ 5 3ao 

, C4=|3aoV^,«o^Y4'7?e 



1^5^26' \/26' 5 V 13^' 

aoV25105T 1747V205 / (±239\/5 + 85^/41) (±13^/5 + 5^41) 1 



75678 



180 



12 



2\/5 
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The corresponding eigenvalues of the scaling matrix at the each fixed points are calculated as 



3 1 

ci; 1 -, - ) , C2; 

55 41 



3 1 



5 -1 



2' 2/ ' V3' 3 



-4 + V63T -4-V63T 



C4; (-5,-11), 



(2.63996417 • • • , -3.63996417 • • •)■ 



Namely, they are classified into the twice unstable fixed points Ci, C2, C5, the once unstable fixed points C3, c± and the 
stable fixed point C4. The critical exponents determined by the once unstable fixed points C3 and c± are respectively 
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5' (-4 + V63T)' 



(68) 



which must be observed most likely in this model. In order to detect other exponents, fine tuning of the initial coupling 
parameters to the twice unstable fixed points is necessary as in the multi-critical behavior of ordinary second-order 
phase transitions. 

Logarithmic dependence of the running coupling constants is controlled by the stable fixed point C4. Since the 
largest eigenvalue at the fixed point is less than —1, the sub-leading contribution is l/(logL)^ in this case. 

It is helpful to study the new RGE for eq.(|67|) numerically in order to better understand our results. Fig.^ exhibits 
the seven fixed points Ci , • • • , c-t in polar coordinates 



with -Tr/2 < (fx 37r/2, < 9 < n. 



a = flo {sin 9 sm(j), sin0 cosc/), cos( 



(69) 





2.5 










.C4 .C.3 ^ 




-1 


0.5 


•C- 1 ,C22 


3^ — ^^^^ 



FIG. 6. Seven fixed points ci, • • • , c± 



The curve in Fig.^ represents the equation a ■ Via) = 0, which corresponds to the set of branch points of the RG 
transformation TZr- AH fixed points belong to the region {a | a • V{a) < 0} because they are on incoming straight 
fiow lines in the original RGE. 

Fig.0 shows the vector field {d(j)/dT, dO/dr) given by eq.(|6|) near the fixed points. From this figure, we find that 
most points in this region go to the outside of the region or to the fixed point C4. 



i i i 4 ^ 4 i i 

i i / / y y / i 



> i 



V V ^ ^ 

T V ^ 



FIG. 7. The vector field defined by the polar-coordinate representation of the new RGE for eq.(^^. The black dots stand 
for the fixed points, which correspond to those in Fig.^. 
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The destination of a trajectory toward the outside of this region is a point on the curve a ■ V{a) = 0. This imphes 
that a flow of the original RGE ( |67|) corresponding to this flow on S approaches the origin once and escapes from it. 
Namely, those points belong to a massive phase. On the other hand, a flow that terminates at the stable fixed point 
C4 corresponds to a flow of the original RGE absorbed into the origin along the incoming straight flow line. The set 
of those points moving toward the fixed point C4 is in a massless phase. 

In general, we have two cases corresponding to massive and massless phases described above. However, there are 
exceptional points which lie on the solutions starting at one of the twice unstable fixed points ci, C2, C5 and arriving 
at one of the once unstable fixed points C3, c±. These exceptional solutions correspond to critical surfaces that give 
the phase boundary determined in the original RGE. More precisely, the set of points on the exceptional trajectories 
is the intersection between a phase boundary determined by the original RGE and the sphere S. It should be noted 
that the phase boundary in the coupling space {(51,52153)} forms a conical surface because V{g) is homogeneous. 




FIG. 8. Solutions numerically computed near the exceptional ones 



Fig.^ depicts a flow numerically solved near the critical surfaces. Using the numerical results Figs.|^ and H, let us 
consider divergent behavior of the correlation length. Suppose that an initial value Oq changes toward the phase 
boundary between Ci and C2 such as the dashed line A. As Oq tends to the phase boundary, the flow starting at ag 
passes near the fixed point C3 and spends a long time there. The result eq.(^) indicates that the critical exponent a 
detected in this case is 3/5. Similarly, if ag varies along lines such as B or C, one observes a — 8/ (—4 + V63T) since 
the solution starting at ao goes through a neighborhood of the fixed points c+ or c_. 

Finally, we comment on the logarithmic dependence of the coupling constants in this example. Figj^ implies that a 



general massless flow is controlled by the twice stable flxed point C4. The universal features considered in sec. IV 



valid if the initial value Oq is sufficiently far from the phase boundary. However, as ao approaches the phase boundary 
from a massless side, the solution is gradually affected by the once unstable fixed points. More profound investigation 
will be needed in this case. 



VI. SUMMARY AND DISCUSSION 

In this paper, we first showed an algebraic way of finding the critical exponent a in eq.(|l|), which was so far com- 
puted by integrating RGE explicitly. The procedure is summarized as follows: 

(i) Derive the new RGE defined in eq.(^) from the original RGE (|^). 

(ii) Find straight flow lines in the original RGE, which correspond to fixed points of the new RGE. 

(iii) Compute the scaling matrix at a fixed point on an incoming straight fiow line and diagonalize it. 

(iv) If the scaling matrix has the unique relevant mode, the correlation length indicates singular behavior by one- 
parameter fine tuning and the exponent cr is equal to the inverse of the relevant eigenvalue. If the scaling matrix has 
multiple relevant modes, we can observe multicritical behavior. 

Second, we derived the logarithmic dependence of running coupling constants in a massless phase where all the 
coupling constants are marginally irrelevant. It was found that the coefficient of the leading log term is universal in the 
sense that it does not depend on an initial value of the running coupling constants, which is the same result as in the 
case of a single marginally irrelevant coupling constant. However, coefficients of subleading terms are non-universal. 
They could disturb the universal nature of subleading terms which come from higher-loop corrections to the original 
beta function. 
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We obtain both results by applying the RG transformation (||) to the original RGE (||), which was motivated by 
the recent developments of RG transformations to non-linear differential equations |^-||] . 

It should be noted that our study is focused when we derive the first result on a flow that goes once toward the 
origin and then leaves for a non-perturbative region. In general, the quadratic differential equation (^) with eq.(|^) 
could have a flow qualitatively different from those we considered. For example, we numerically find that the equation 




FIG. 9. Flow in eq.Q with g = 2 



Namely, the equation has solutions that first escape from the origin and then turn back to it. In this case we cannot 
directly apply the result eq.(^). Even in that case, as we show in the following, the beta function in eq. (0) is helpful 
for understanding a qualitative picture of these solutions. In fact, we first notice that there are only four straight 
flow lines on f;i = or on 02 = 0. Next let us compute the scaling matrix in eq.(^) at the fixed points (±ao, 0) and 
(0, ±ao) of the new RGE ([lOl). We find that it has eigenvalues 1 — g at (±ao, 0) and at (0, ±ao). If g > 1, solutions 
for the RGE near (±ao,0) is absorbed into the fixed points. The integral curve does not change the direction at 
the other fixed points (0, iag) because the eigenvalue of the scaling matrix vanishes at those points. Therefore, by 
continuity, we conclude that there must be at least two branch points on S. 




FIG. 10. Flow in the new RGE for eq.([7o|) with q = 2. The white circles correspond to turning points. 



Moreover, the solution must escape from the branch points. As we pointed out in the last part of II A, it means 
that eq.(|70|) has a flow that first leaves from the origin and turns back to it, as depicted in Fig. ||. 

The study of the universality classification of infinite-order phase transitions is now in progress. In two dimensions, 
it should contribute to the analysis of the c > 1 CFT. In higher than two dimensions, a nontrivial exponent in an 
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infinite-order phase transition might be observed experimentally in some phenomena. 
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APPENDIX A: RELATIONSHIP BETWEEN THE EIGENVALUES OF A{a*) AND B{a*) 



In this appendix we prove eq.(36). The claim is that a set of the eigenvalues of B{a*) in cartesian coordinates is 
equal to that of A{a*) in polar coordinates plus extra zero eigenvalue. 

To show that, we add a/ a = e„ to the basis {eQ}i<Q,<„_i of the tangent space at a(r) g S. Then the set {ei}i<i<„ 
becomes an orthonormal basis of the n-dimensional space of coupling constants. Consider the n x n matrix 



Tij{a) = f.j{a){ei,ej), 



(Al) 



where {ei}i<i<„ is the orthonormal basis that defines the cartesian coordinates ((^i, • • • , (7„) and the bracket {x,y) 
means the inner product. The function fj{a) is given in eq.jl^) for 1 < j < n — 1. In addition, 



fnia) = 



da 



da 



= 1. 



(A2) 



Since (e^, ej) forms an orthogonal matrix, wc immediately have the inverse of T as 

T-,\a) = Ma)-'ie.„ek). 
Now we examine the form of the n x n matrix T^^BT. As the first step, let us compute T^^B: 

fe=i fe=i °' 

Here we have used l3{a*) = in the second equality. Next, according to eq.(p^), we find that can be written as 



(A3) 



(A4) 



„ , da , dai , 
Tij^ieug^)^^ (l<J<n-l). 



d9i 



Using eqs.( |A4|) and (|A5|) , we get 



Y: Tr^\a*)Bu{a*)TiAal ^ E = 

l<k,l<n l<k,l<n ' ' 

for 1 < i < n and 1 < j < n — 1. 

The result indicates that T^^BT has the following form: 



(A5) 



(A6) 



(T-^BT) (a*) 



A{a*) 



(A7) 



which proves eq.(M). Note that the last row vanishes because /3n(a) = (/3(a), e„) = for all a. 
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